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INTRODUCTION 



DATASET AND APPARATUS 



One decade after the first observation of exclusive 
B K*j decays by CLEO in 1993 Q, this process con- 
tinues to be a subject of considerable interest. The size 
of the decay rate itself provides only a mild constraint 
on extensions to the Standard Model (SM), because SM 
predictions for exclusive rates suffer from large (^--^30%) 
and model dependent form factor uncertainties 0. Of 
more interest are asymmetries, where theoretical uncer- 
tainties largely cancel. The isospin asymmetry between 
the charged and neutral B K*j decay widths is pre- 
dicted to be +5 to 10% in the SM, while in some SM 
extensions it may have an opposite sign j^]. A measure- 
ment of the partial rate asymmetry between CP conju- 
gate modes is another interesting subject; here the SM 
expectation is much less than 1% and any large asym- 
metry would be an indication of non-SM effects. In this 
report, we present new measurements of the B K*j 
branching fractions, and isospin and charge asymmetries. 



The data sample used in this analysis contains (85.0 ± 
0.5) x 10^ B meson pairs, corresponding to an integrated 
luminosity of 78 fb^^, collected at the T(4S') resonance 
by the Belle detector at the KEKB storage ring. KEKB 
is a double-ring asymmetric-energy e'^e~ storage ring 
(3.5 GeV on 8 GeV) Q. We also use an off-resonance 
data sample of 8.3 fb~^ collected at a center-of-mass 
(CM) energy that is 60 MeV below the T{AS) resonance. 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a three-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel- like arrangement of time-of- flight scintil- 
lation counters (TOP), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside 
a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to identify 
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muons fKLM) . The detector is described in detail else- 
where Isl. 



SIGNAL RECONSTRUCTION 

The analysis is performed by reconstructing B meson 
candidates that include a high energy primary photon 
and a K* resonance reconstructed in one of four final 
states: K+Tr~, K°Tr", K°tt+, and K+ir^. Here and 
throughout this report, K* denotes the if* (892), and 
the inclusion of charge conjugate modes is implied unless 
otherwise stated. 

Photon (7) candidates are reconstructed from isolated 
clusters in the ECL that have no corresponding charged 
track, and a shower shape that is consistent with that 
of a photon. The photon energy is calculated from the 
sum of energies in crystals with more than 0.5 MeV en- 
ergy deposited around the central cell. Photons in the 
energy range 1.8 GeV < E* < 3.4 GeV in the barrel 
region of the ECL (33° < 9j < 128°) are selected as pri- 
mary photon candidates from B decay; here E* is the 
photon energy in the CM frame and O-y is the polar an- 
gle in the laboratory frame. (We use variables calculated 
both in the CM frame and laboratory frame: variables 
defined in the CM frame are labeled with an asterisk.) In 
order to reduce the backgrounds from tt" and r] mesons 
from continuum light quark-pair production (e^e~ — > qq, 
q = u,d, s, c), we impose two additional requirements on 
the primary photon. One is the explicit removal of tt" 
(77) candidates by requiring the invariant mass of the 
primary photon and any other photon with an energy 
greater than 30 MeV (200 MeV) to be outside of a win- 
dow of ±18 MeV/c^ (±32 MeV/c^) around the nominal 
TT*^ (77) mass. These correspond to ±3cr windows, where 
a is the mass resolution. This set of criteria is referred 
to as the n'^/rj veto. The other requirement is the re- 
moval of the clusters that are not fully consistent with 
an isolated electromagnetic shower. We require the ratio 
of the energy deposition in 3 x 3 cells to that in 5 x 5 
cells around the maximum energy ECL cell of the cluster 
(^9/25) to be greater than 0.95, which retains 95% of the 
signal photons. 

Charged pions (7r+) and kaons {K~^) are reconstructed 
as tracks in the CDC and SVD. The tracks are required 
to originate from the interaction region by requiring that 
they have radial impact parameters relative to the run- 
averaged measured interaction point of less than 1.5 cm. 
We determine the pion (L^) and kaon (Lk) likelihoods 
from the ACC response, the specific ionization (dE/dx) 
measurement in the CDC and the TOF flight-time mea- 
surement for each track, and form a likelihood ratio 
Lk/tt — Lk/{Lt^ + Lk) to separate pions and kaons. We 
require Lxjix > 0.6 for kaons, which gives an efficiency of 
86% for kaons, and L^/tt < 0.9 for pions, which gives an 
efficiency of 96% for pions. In addition, we remove kaon 



and pion candidates if they are consistent with being elec- 
trons based on the ECL, dE/dx, and ACC information. 

Neutral pions (tt^) are formed from two photons with 
invariant masses within ±16 MeV {3a) of the 7r° mass; 
the photon momenta are then recalculated with a ir^ mass 
constraint. The tt" mass resolution is better than that 
for the TT^'/yy veto where the photon energies are highly 
asymmetric. We require each photon energy to be greater 
than 50 MeV, and the cosine of the angle between the 
two photons (cos6'.y^) to be greater than 0.5. This angle 
requirement is almost equivalent to selecting tt'^s with 
momentum above 0.5 GeV/c; it retains about 90% of the 
signal 7r°s while rejecting 43% of the 7r° candidates in the 
background. 

Neutral kaons (Kg) are reconstructed from two oppo- 
sitely charged pions that have invariant masses within 
±10 MeV (3cr) of the Kg mass; the pion momenta are 
then recalculated with a Kg vertex constraint. We im- 
pose additional criteria based on the radial impact pa- 
rameters of the pions {6r), the distance between the 
closest approaches of the pions along the beam direc- 
tion (Sz), the distance of the vertex from the interaction 
point (l), and the azimuthal angle difference between the 
vertex direction and the Kg momentum direction ((50). 
These variables are combined as follows: for p{Kg) < 
0.5 GeV/c, Sz < 8 mm, Sr > 0.5 mm, and S(j) < 0.3 rad 
are required; for 0.5 GcV/c < piK*^) < 1.5 GeV/c, Sz < 
18 mm, Sr > 0.3 mm, Scj) < 0.1 rad, and I > 0.8 mm 
are required; and for p{Kg) > 1.5 GeV/c, Sz < 24 mm, 
Sr > 0.2 mm, S(p < 0.03 rad, and / > 2.2 mm are re- 
quired. This set of criteria retains about 80% of the sig- 
nal K^. 

We form a B candidate from a primary photon can- 
didate and a K* candidate, which is a Kn system with 
an invariant mass M{K'k) within ±75 MeV/c^ of the K* 
mass. In order to separate the B candidate from back- 
grounds, we form two kinematic variables: the beam- 
energy constrained mass M^c = \/(-E'bcam/c^)^ ~ IPb/^P 
and the energy difference AE — Eg — E^^^^^, where 
-^boam the beam energy, and Eg and pg are the 
energy and momentum, respectively, of the B candi- 
date in the CM frame. The energy Eg is calculated as 
Eg = E* + Ej^,; the momentum pg is calculated with- 
out using the absolute value of the photon momentum 
according to 

Pb = PK' + X (^bcam ^E^'), (1) 

since the K* momentum and the beam energy are de- 
termined with substantially better precision than that of 
the primary photon. 

We use Mbc as the primary distribution to extract the 
signal yield. For modes without a tt*^, we use a Gaussian 
function with a width of (2.73 ± 0.04) MeV/c^ to model 
the signal; for modes with a tt'^, we use an empirical for- 
mula to reproduce the asymmetric ECL energy response 
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(known as the Crystal Ball line shape (g), whose effec- 
tive width is (3.35 ± 0.10) MeV/c^. The peak positions 
and widths are primarily determined using Monte Carlo 
(MC) samples, and corrected for the measured differences 
in the beam-energy and its spread between data and MC 
using a i?" D^tt" sample. The AE signal distribution 
also has a large tail on the negative AE side due to the 
asymmetric ECL energy response. For the modes with- 
out a 7r°, we use a Crystal Ball line shape; for the modes 
with a 7r°, we convolve an additional Gaussian resolution 
function to describe a broader width and add a broad 
Gaussian component for the small tail in the positive AE 
side. These shapes are determined using MC samples. 
We select candidates with -200 MeV < AE < 100 MeV 
to accommodate the asymmetric AE signal shape. We 
define a AE sideband as 100 MeV < A^; < 400 MeV, 
where no signal is expected, to study the Mbc distri- 
bution of the background. There is no background 
that makes a peak in this AE sideband. We require 
Mbc > 5.270 GeV/c^ when the AE distribution is exam- 
ined. We define 5.227 GeV/c^ < ilfbc < 5.263 GeV/c^ as 
an Mbc sideband that is used to study the AE distribu- 
tion of the background. 



BACKGROUND SUPPRESSION 

The main background source is continuum qq produc- 
tion including the initial state radiation process e'^e" —> 
qqj. We reduce this background by exploiting the topo- 
logical event shape differences: B meson pairs decay al- 
most at rest in the CM frame and thus the final state 
particles are distributed nearly isotropically; qq pairs are 
produced back-to-back with multi-GeV/c momenta in 
both hemispheres and, thus, tend to be more two-jet like. 

We define a Fisher descriminant (F) from modified 
Fox- Wolfram moments la, 



(2) 



(=1 



where a;. Pi are coefficients that are selected to provide 
the maximum discrimination between the signal and the 
continuum background. The modified Fox- Wolfram mo- 
ments are defined as 



(3) 



J2^.J \P^\ 



where the indices i, j indicate the charged tracks (with 
a 7r+ mass hypothesis) and photons that are not used to 
form the B candidate, and the index 7 corresponds to 
the primary photon. The variables 9^^ and 9'^^ are the 



opening angles between two momentum vectors, and Pi 
is the Z-th Legendre polynomial function. The momenta 
(p ') and angles {9') are calculated in the candidate B rest 
frame (denoted as primed variables), since the selection 
efficiency with this variable has a smaller correlation with 
Mbc than that calculated in the CM frame. 

As an additional discriminant, we use the cosine of 
the CM polar angle of the candidate B flight direction, 
cos 9g . The cos 9^ distribution is 1 — cos^ 9g for B pro- 
duction from e+e" T(4S'), and is found to be flat for 
the continuum background. 

We combine these two discriminants into a likelihood 
ratio, 



■^cont — 



Ls + Lb ' 
Lb 



(4) 



pF' V pcosB 



where and P"^"^^ are the probability density functions 
(PDF) for the Fisher and the B flight direction, and the 
indices S and B denote the signal and background. For 



we use |(1 — cos"^ 9g) and ^, respec- 



P|°"^ and Plj"*^^ 
tively. For Pg and Pg , we model the shape by fitting 
the signal and continuum background MC distributions 
with asymmetric Gaussian functions for each of the four 
B K*j channels. 

The value of £cont ranges between and 1. We op- 
timize the minimum £cont requirement to provide the 
largest value of Ns/\^Ns + Nb, where Ns and Nb are 
the expected signal and background yields for 85.0 x 10^ 
B meson pairs assuming previo usly measured B K*j 
branching fraction values |^, Hol lllj . Although there is a 
slight mode dependence in the optimal value, we apply 
the same requirement, £cont > 0.65, which is close to the 
optimal point for each mode. This requirement retains 
73% of signal events while rejecting 90% of continuum 
background events. 

The remaining continuum background is distinguished 
by fits to the Mbc distribution. The continuum back- 
ground is modeled with a threshold function (known as 
the ARGUS function [12), 



/cont(Mbc) 



N X Mbc X 



X exp < a 



E, 
Mbc 

"^bcam 



be 



beam 
2 



(5) 



where A'^ is a normalization factor and a is an empiri- 
cal shape parameter. We determine the shape param- 
eter from the AE sideband data, since there is no sig- 
nificant difference between the background shapes in the 
AE sideband and the signal region for off-resonance data 
and MC events. The results are consistent for data and 
MC, and for the different B ^ K*j channels. We use the 
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same background shape parameter a for all the B K*^ 
channels. For AE', the background shape is parameter- 
ized as a linear function and determined from Mbc side- 
band data. 

Major background contributions from B decays are 
from cross-feeds between charged and neutral B K*j 
decays, B {K*tt + Kp)-f ^ K*rj Q and the 

unmeasured mode B K*it^ |15| for which we assume 
half of the upper limit as the branching fraction with a 
100% error. These backgrounds peak in Mbc around the 
signal with a slightly larger tail, and have a broad peak in 
/S.E at negative values. We model these backgrounds and 
other B decay backgrounds with a smoothed histogram 
generated from a large MC sample. 

SIGNAL EXTRACTION 

We extract signal yields in each of the four final states 
using a one-dimensional binned likelihood fit. The Mbc 
distributions are modeled as a sum of three components: 
the signal, the continuum background and the B decay 
backgrounds that are described in the previous sections. 
Figure ^ shows the result of the fits; clear signals are 
seen in all four final states. The size of the B decay back- 
ground component, which is seen as a slight enhancement 
of the background shape under the signal peak, is fixed in 
the fitting procedure. We vary the B decay background 
components by the errors on their branching fractions to 
evaluate the systematic error due to their uncertainties. 
We also vary the continuum background shape, i^bcamj 
the B meson mass, and the Mbc resolution by their errors 
to evaluate the systematic errors on the signal and back- 
ground PDFs. We use a quadratic sum of the variations 
in the signal yield with these tests as the systematic error 
on the signal yield. 

A similar fitting procedure is performed for the A_B dis- 
tributions as shown in Fig.|21 The IS.E yields are obtained 
by integrating the fit results from —0.2 GeV to 0.1 GeV 
to allow a comparison with the Mbc results. The fit re- 
sult for each mode is in agreement with the Mbc result 
as given in Table HI In this analysis, the Ai? distribu- 
tion is not as reliable as Mbc since we have to consider a 
wide \E range where the background contribution from 
B decays is larger than that for Mbc, and /S.E shapes for 
both signal and background may have large uncertainties 
that are not fully evaluated due to lack of suitable control 
samples. Thus we base our signal yields on the Mbc fits. 

The M{Kit) invariant mass spectrum, before apply- 
ing the \M{K'k) — Mk* \ < 75 MeV/c^ requirement, gives 
discrimination of the K* signal from resonances such as 
i4r|(1430) and iir*(1410), or a non-resonant Knj compo- 
nent. In order to examine the spectrum, we divide the 
data below M{KTr) = 2.0 GeV/c^ into 50 MeV/c^ wide 
bins and extract the signal yield for each bin from a fit 
to the Mbc distribution, using the same fitting procedure 



described above. We veto the K 7r+ and D° 

K^TT^ contributions from W D^n" and W D^rj 
backgrounds by requiring \M{Kt:) — Mdo \ > 20 MeV/c^ 
for the iC*°7 modes. Other B decay backgrounds that 
may peak in Mbc are included as a background compo- 
nent in each fit. A fit is then performed to the M{Ktt) 
spectrum using a sum of K* , i5r|(1430) and K*{\A\Q) 
resonances and a non-resonant component without tak- 
ing into account possible interference effects. The fit 
results are shown in Fig. |3 For these resonances, we 
use relativistic Breit-Wigner functions with the nominal 
masses and widths convolved with a Gaussian resolution 
function. In addition to the K* peak, we find a signif- 
icant X2(1430) peak, while the i^*(1410) component is 
consistent with zero. The K* signal yields within the 
±75 MeV/c^ window are consistent with the Aibc results 
as given in Table In all cases, the background con- 
tributions from K*{141Q) and i4:2*(1430) within the K* 
mass window are less than one event including their er- 
rors. The non-resonant contribution is modeled with the 
inclusive Xg mass spectrum of Ref. We find 5.9 ±2.3 
and 6.4± 1.9 events under the peak of K*^ and if*"'", cor- 
responding to (1.2 ±0.5)% and (2.4 ±0.4)% of the signal 
yields, respectively. We include these yields into the sys- 
tematic errors instead of subtracting the contributions, 
since they could also be due to a bias from other high 
mass resonances such as i(r*(1680) that hardly contribute 
to the K* mass peak. 

The K* decay helicity angle, 6'hci, defined as the an- 
gle between the kaon and the B meson directions in the 
rest frame of the Kt: system, provides discrimination of 
the spin-1 signal from other spin states. The cos^hci dis- 
tribution, shown in Fig. ^ is obtained by dividing the 
data into bins of 0.1 in cos^hci and fitting their Mbc dis- 
tributions. The B decay backgrounds that peak in Mbc 
and follow a cos^ 0hci structure due to the pseudoscalar 
to pseudoscalar-vector nature of the decay are included 
as a background component in each fit. The K*^ signal 
has a 1 — cos^ 0hci distribution. In order to model a slight 
distortion due to non-uniform tracking, particle identi- 
fication and TT^ /Kg reconstruction efficiencies, we use a 
fourth order polynomial function that is constrained to 
zero at cos ^hei = ±1. We also add a flat component mod- 
ified with the same slight distortion, which turns out to 
be consistent with zero. The fit result is in agreement 
with the spin-1 signal. 

We obtain a consistent set of signal yields from three 
different distributions, Mbc, Ai?, M{Ktt) as summarized 
in Table m Both the M{Ktt) and cos6'hoi distributions 
suggest that the contributions from other resonances and 
non-resonant decays can be neglected. We conclude that 
the signal yields obtained with the Mbc fit are essentially 
entirely due to B ^ K*^ decays. 
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BRANCHING FRACTIONS 

The reconstruction efficiencies are primarily obtained 
from signal MC samples. The selection criteria are di- 
vided into ten categories, and the systematic error for 
each of them is evaluated with an independent control 
sample. The results are summarized in Table Hll 

The uncertainty in the photon detection efficiency is 
evaluated with a sample of radiative Bhabha events. For 
the tracking efficiency, we quote an error from a com- 
parison of the partially reconstructed D*'^ D^tt^, 
KgTT~^Tr~ yield with the fully reconstructed one. 
For the charged kaon identification, we evaluate the sys- 
tematic error from a comparison between the efficien- 
cies obtained from kinematically selected D'^tt^, 
D° K^Tr+ and cj) K+K^ decays in data and MC. 
Similarly for the charged pion identification, we com- 
pare the efficiency for the same D*+ sample with that 
obtained from Kg tt+tt^ to evaluate the system- 
atic error. The uncertainty in the Kg reconstruction 
is obtained from a comparison between _D+ -^s""^ 
with D+ -> K'Tr+TT+ in D*+ D+tt° decays. The 
uncertainty in the tt" reconstruction efficiency is evalu- 
ated from a comparison of 77 — > Tr^Tr'^Tr" to 77 — > 77 and 
77 7r+7r^7r° in data and MC. The efficiencies for the 
77*^/77 veto and the likelihood ratio requirement are eval- 
uated together, using a B Dtt^ control sample that 
includes the decay channels D'^tt^ , K^tt^ 

and ^ D+TT", D+ K-tx+h+. We calculate Mbc as 
we do for B K*'f, i.e. without using the absolute value 
of the 7r~ momentum. For the Tr^/Ty veto, we assume the 
7r~ is a massless particle and scale the momentum by a 
factor of 1.1 to make the average tt~ momentum equal 
that of the photons for B K*^. We then combine 
this massless 7r~ with all photon candidates and reject 
the event if the same tt^/t? veto criteria are satisfied, and 
compare the results between data and MC control sam- 
ples. The effect of the M{Ktt) requirement is evaluated 
by taking into account the K* form factor and the de- 
tector resolution effects into the Breit-Wigner shape. 

Using these reconstruction efficiencies, we obtain the 
branching fractions for each of the four modes as sum- 
marized in Table ITTll We add the signal yields and the 
efficiencies for two modes of each of neutral and charged 
B decays, and obtain 



B{B 
B{B 



K*°"f) =(4.01 ± 0.21 ± 0.17) X 10-5 
K*+j)^Ia.25 ± 0.31 ± 0.24) x 10^^^ 



(6) 



where the first and the second errors are statistical and 
systematic, respectively. We assume an equal production 
rate for B°B^ and B^ B~ from the T(4S') resonance. 
The isospin asymmetry. 
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(Tg+/rgo)g(j3" -> JC*°7) - B{B+ ^ K*+^) 
{tb+/tb»)B{B° K*Oj)+B{B+ K*+-iy 



is then calculated from these results. We use the world 
average value of tb+ /tbo = 1.086±0.017 0|. We assume 
that the systematic error on the photon detection cancels. 
The systematic error on the £cont and tt^/t^ veto require- 
ments is estimated to be 0.013 in Ao+ from a comparison 
of the B^ and B^ subsets of the B Dtt^ control sam- 
ple. We find the systematic error due to the M{Ktt) 
requirement is negligible. Correlations between the sys- 
tematic errors for the charged pion and kaon tracking and 
particle identification, and the tt" and Kg reconstruction 
efficiencies are taken into account. Systematic errors on 
the fitting procedures are assumed to be uncorrelated. 
The result is 

Ao+ = +0.012 ± 0.044(stat) ± 0.026(syst), (8) 

which is consistent both with the SM prediction and 
no asymmetry. Although the systematic errors in the 
branching fractions are almost as large as the statistical 
errors, the systematic errors largely cancel in Ao+. 

If we allow the 5+ to B° production ratio {f+/fo) to 
deviate from unity, the value of Ao+ shifts approximately 
by ^(U/fo - !)■ The value of /+//o = 1.044 ± 0.050 in 
Ref. m gives Ao+ = +0.034±0.044(stat)±0.026(syst)± 
0.025(/+//o). The conclusion above is therefore un- 
changed, although the result is shifted closer to the SM 
prediction. 



(7) 



SEARCH FOR PARTIAL RATE ASYMMETRY 

We define the partial rate asymmetry between CP con- 
jugate modes (except for the Kgir'^j mode) as 

r(B^K*-f)-T{B^K*j) 

Anp — — = = 

V(B K*-f) + r{B K*-f) 

1 N(B -> K*-f) - N{B K*-f) 
1 - 2w ^ N(B -> K*-f) + N{B K*-f) ' 

(9) 

where N is the signal yield, w is the wrong-tag fraction, B 
indicates either B° or K* indicates K*°{-^ K+n^) 
or K*+{-^ iCg7r+,X+7r°), and B, K* are their conju- 
gates, respectively. 

The wrong-tag fraction is negligible for K*'^, and is 
0.9% for K*'^ due to the double mis-identification of 7r+ 
as K'^ and K~ as 7t~ . The wrong-tag fraction is obtained 
from the signal MC; we neglect the small error on this 
fraction. 

Possible detector and reconstruction biases are studied 
with an inclusive K* sample. We compare the yield of 
K* and K* and find no significant difference in Ak* = 

ii-2w)N{?']lmK'] - conclude there is no bias and 
assign systematic errors of 0.007 for K*'^ A'+tt^ and 
K*+ K^n^, and 0.015 for K*+ JsT+tt". The sys- 
tematic error on the £cont and 7r'^/77 veto requirements is 
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estimated to be 0.007 in Acp from a comparison of the 
B and B subsets of the B D-k~ control sample. 

We divide the data shown in each of Fig. H^a-c) into 
CP conjugate modes. We fit the Mbc distributions for 
the six modes separately, using the same fitting proce- 
dure used in the branching fraction measurement. The 
fit results are shown in Fig. [31 and summarized in Ta- 
ble |W] The errors on the yield extraction are assumed 
to be uncorrelated. By assuming the partial rate asym- 
metry is equal for charged and neutral B decays, we add 
the signal yields for B and B with a small correction due 
to the wrong-tag fraction, and obtain 

Acp{B K*"f) = -0.015 ± 0.044(stat) ± 0.012(syst). 

(10) 

Here the systematic error includes the errors on Apc* , 
£cont and the yield extraction. 

CONCLUSIONS 

We have presented branching fraction, isospin asym- 
metry and CP asymmetry measurements for the radia- 
tive decay B K*^ using 85.0 x 10^ B meson pairs. 
The branching fraction results are consistent with pre- 
vious Belle [3 results, and also with the CLEO ^3 
and BaBar |ll| results, with the errors improved by 
a factor of two. The Ktt mass spectra and the de- 
cay helicity angle distribution in the K* mass region 
{\M{K'k) - MkA < 75 MeV/c^) are consistent with the 
dominance of B — > K*"f without other contributions such 
as 1410)7 or non-resonant Knj decays. We measure 
an isospin asymmetry which is consistent with zero; with 
the current precision, our result agrees with both the SM 
prediction and new physics scenarios which have the op- 
posite sign. For the partial rate asymmetry between CP 
conjugate modes, we obtain a result which is also con- 
sistent with zero. For both of these asymmetries, the 
systematic errors are much smaller than the statistical 
errors, and hence we can expect further improvements 
with larger data samples. 
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FIG. 1: Fit results for the beam-energy constrained mass dis- 
tribution for the (a) K+tt''), (b) Kg-K+^i, (c) K+tt^-), and 
(d) T^gTT^ modes. The sum of the signal and the background 
components are shown in the solid curves, while the dotted 
and dashed curves represent the total backgrounds and the B 
decay backgrounds, respectively. 
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FIG. 2: Fit results for the AiJ distributions. The sum of the 
signal and the background components are shown in the solid 
curves, while the dotted and dashed curves represent the total 
backgrounds and the B decay backgrounds, respectively. 



TABLE I: Summary of the signal yields. 

Mode Mbc Ag M(K-k) 

K+'R-'y 450.1 ± 24.2 ± 6.1 453.2 ± 27.4 

-K'g7r°7 23.8 ±6.4 ± 1.0 23.2 ± 6.2 

K%Ti+'y 145.0 ± 13.7 ±3.6 147.8 ± 15.4 

_K'+7r°7 129.1 ±14.7 ±5.5 120.0 ± 14.8 

if*°7 473.9 ± 25.0 ± 6.2 476.4 ± 28.1 480.6 ± 25.7 
K'+'y 274.1 ± 20.1 ± 6.6 267.8 ± 21.4 270.9 ± 19.9 
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FIG. 3: Fit results for the Kit invariant mass distributions for 
the (a) sum of K'^-K~'y and K'g-K^'^, and (b) sum of K%-k^^ 
and K channels. The sum of the signal and the back- 
ground components are shown in the solid curves, while the 
dotted curves represent the background components. 
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FIG. 4: Fit results for the cosine of the helicity angle (cos Shci) 
distribution for the sum of all four B K*y channels. The 
solid line shows the sum of the signal and the flat component; 
the ±lo" bounds of the latter component are shown as the 
dotted lines. 



TABLE 11: Reconstruction efficiencies and their systematic uncertainties. 
K+TT--/ Ji'g7r+7 K+TT°'y A'g7r°7 



Reconstruction efficiency (12.83 ± 0.54)% (3.96 ± 0.25)% (3.63 ± 0.19)% (1.09 ± 0.07)% 
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FIG. 5: Fit results for the beam-energy constrained mass dis- 
tributions for the search for partial rate asymmetry. The sum 
of the signal and the background components are shown in the 
solid curves, while the dotted and dashed curves represent the 
total backgrounds and the B decay backgrounds, respectively. 



TABLE III; Results for the signal yields, efficiencies and 
branching fractions (B). 



signal yield efficiency B (xlO 
K*°^ 473.9 ± 25.0 ± 6.2 13.92 ± 0.58 4.01 ± 0.21 ± 0.17 
K*+y 274.1 ±20.1 ±6.6 7.59 ±0.39 4.25 ± 0.31 ± 0.24 



(A"+7r~)7 450.1 ± 24.2 ± 6.1 
{K^Tv^h 23.8 ±6.4 ±1.0 
(/i'g7r+)7 145.0 ± 13.7 ±3.6 
{K+TT°)'y 129.1 ±14.7 ±5.5 



12.83 ± 0.54 4.13 ± 0.22 ± 0.18 
1.09 ±0.07 2.57 ±0.69 ±0.20 
3.96 ±0.25 4.31 ±0.41 ±0.29 
3.63 ±0.19 4.19 ±0.48 ±0.28 



TABLE IV: Results of the partial rate asymmetry search. 

K* mode NjB K'-y) N{B K*'y) Acp 

K*° K+TT- 218.5 ± 16.8 ± 3.0 231.6 ± 17.4 ± 3.0 -0.030 ± 0.055 ± 0.014 
K*+ ^ K°s7v+ 79.2 ± 10.0 ± 1.9 65.8 ±9.4 ±1.9 ±0.094 ± 0.094 ± 0.021 
K*+~*K+Tv'' 68.7 ±11.6 ± 2.2 80.1 ± 12.4 ± 2.5 -0.078 ± 0.113 ± 0.028 
Combined {K*+) ±0.007 ± 0.074 ± 0.017 



Combined (all) 



-0.015 ±0.044 ±0.012 



